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Abstract 

We examine the vacuum overlap formula for the two-dimensional SU(2) Wess- 
Zumino term in lattice regularization. Perturbatively it reproduces the Wess-Zumino 
term correctly in the continuum limit and yields the IR fixed point in the beta function 
of the chiral model. Nonperturbatively it shows a sharp Gaussian distribution for the 
SU(2) chiral field configurations in the scaling region, where smooth configurations 
dominate even in the symmetric phase due to asymptotic freedom. Crossover is sharp 
from the strong coupling region where the Wess-Zumino term fluctuates hard and the 
species doublers' contribution is suspected to affect it. 
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§1. Introduction 

The vacuum overlap formula§)>i>>i'i)>i>'§ gives a well-defined method to calculate the 
chiral fermion determinant in lattice regularization. By this formula, one first considers two 
Wilson fermion fields which posses positive mass and negative mass, respectively. They are 
assumed to couple to the link variables of given background gauge field. We next solve the 
ground states of the two fields. Then the overlap of the two ground states gives the chiral 
determinant. There remains an ambiguity in the complex phase of the determinant. It is 
fixed by the convention that the inner product of each ground state with the corresponding 
free ground state is real and positive. 

The complex phase of the vacuum overlap is able to reproduce the full complex phase 
of chiral determinant correctly for a smooth background gauge field. The complex phase is 
gauge noninvariant and its variation under a gauge transformation was shown to give the 
consistent anomaly. In the continuum U(l) chiral gauge theory which is defined 
on the two-dimensional torus, the chiral determinant with the uniform background gauge 
potential (toron) can be calculated exactly. The vacuum overlap reproduces this exact 
result. § 

Actually, this lattice definition of the complex phase of chiral determinant is closely 
related to the i] invariant definition in the continuum theory. 0) The vacuum overlap formula 
can be derived from the (2n + 1) dimensional system of Wilson fermion with kink mass0) by 
taking infinite limit of the extra dimension, provided that appropriate Pauli-Villars (Wilson) 
fields are included. ) >lll)>ll ) For this limiting procedure to be well-defined, it also needs that 
the background gauge field is topologically trivial. B* By virtue of the Pauli-Villars fields, 
we may take the naive continuum limit a — > in this formula keeping the kink and Pauli- 
Villars masses finite. Then it reduces to the formula which defines the rj invariant. The only 
technical difference is that the smooth interpolation function of two gauge fields is replaced 
by the sharp step function at the boundary of the wave guide of gauge interaction. 

For generic gauge fields on the lattice, however, several authors E1)>E1)>0>'!)>@) argued that 
the nature of gauge degrees of freedom may spoil the chiral structure of the formula. For 
the (2n + 1) dimensional Wilson fermion with the kink mass and the wave guide of gauge 
interaction, 0) it was shown© that the large fluctuation of gauge freedom induces another 
light fermion mode of opposite chirality at the boundary of the wave guide. Such a light 
mode could also give the contribution to the complex phase of chiral determinant. 

Singular gauge transformations can also affect the complex phase of the vacuum overlap. 
In the continuum U(l) gauge theory which is defined on the two-dimensional torus, we 
can consider the singular gauge transformation which wipe out the toron into the singular 
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potential with delta functions. Then, it was showncP that the gauge invariance breaks by the 
singular gauge transformation even in the anomaly free U(l) chiral gauge theory, provided 
that the background toron is large and the anti-periodic boundary condition is chosen for all 
fermions. The vacuum overlap reproduces the result as well. The authors of ref.H* claimed 
that it is these singular gauge transformations which spoils the chiral phase when averaged 
on the gauge orbit. 

Therefore, it is a crucial problem in this approach to control and reduce the fluctuation 
of gauge degree of freedom. In this paper, we examine the gauge freedom of the two- 
dimensional SU(2) gauge field. In order to control and reduce the fluctuation, we introduce 
the kinetic term for the variable of gauge freedom, which consists of the nearest neighbor 
coupling. We calculate the Wess-Zumino term by the vacuum overlap formula and measure 
its distributions. The whole system we consider here is the lattice counterpart of the SU(2) 
Wess-Zumino- Witten model off the critical point. 

This type of kinetic term for gauge freedom was examined in the U(l) case.© It was 
shown that the reduction of the fluctuation is not sufficient near the critical point of the 
corresponding X-Y model. On the contrary, in the case of SU(2) chiral field considered here, 
the coupling in the critical region is weak because of the asymptotic freedom. It is expected 
that the fluctuation can be reduced in this region even thought the SU(2) chiral symmetry 
is realized linearly. 

This paper is organized as follows. In section 2, the lattice counterpart of the two- 
dimensional SU(N) Wess-Zumino- Witten model is examined in the lattice perturbation the- 
ory. We show that the vacuum overlap formula reproduces the Wess-Zumino term correctly 
in the continuum limit. We also calculate the beta function of the model and show that 
it has the nontrivial IR fixed point, as in the continuum theory. In section 3, the Wess- 
Zumino term is examined numerically. The SU(2) chiral field is generated by the Monte 
Carlo method (Cluster algorithm©) with the action of the nearest neighbor coupling for 
several values of the coupling constant. With these configurations, we calculate the complex 
phase of the vacuum overlap and its distribution is obtained. In section 4, we give a discus- 
sion concerning the numerical estimate of observables in the Wess-Zumino- Witten model by 
incorporating the imaginary Wess-Zumino action. 



*) For the small toron, thexhiral determinant in the continuum limit is gauge invariant even under the 
singular gauge transformation, u Choosing such small toron, the gauge integration of chiral determinant was 
performed. EiP The result shows that the large gauge fluctuation affects the chiral phase even for the small 
toron and seems to cause it to vanish. 
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§2. Two-dimensional Wess-Zumino term on the Lattice 

2.1. Vacuum overlap formula and Wess-Zumino term 

The definition of the vacuum overlap formula for the lattice regularized chiral determinant 
in In dimension is given byS) 



det(C) rcg = , 7 , , , V (2.1) 

\\\—\—)v\ x\— 

|±)[/ are the vacuum state of the second quantized Hamiltonian H ± , which defined by 



H\U) 




(2.2) 




ram rra/3j 




H±(U) 


= 75[Tm + C + S], 


(2.3) 




= 7^ (Tm) a/3 {fim,n+lJ,Un,U ^ra,m+/i£^rra,/j) 5 


(2.4) 


Bnai,mf3j 


7} ^ y <^a/3(2<$ra,m ^m,n+fiUn,fjb ^ra,m+A^m,/j) ) 


(2.5) 



where we set the lattice spacing a to unity, mo is a bare mass parameter , C and B are 
the massless Dirac operator and the Wilson term respectively. a na i, a nai are field operators 
satisfying the following anti-commutation relations : 

The vacuum overlaps, u(~\+)u etc., can be expressed by the one particle wave functions 
^(n, U), which are the eigenfunctions of H ± {U) with positive eigenvalues labeled K: 



v 



'-\+)u = det (£V£(n, tf)ty£(ra,E0) , (2-7) 
,(±|±)! = det fev£(n,tf)ty|,(M) ] , (2-8) 



where color and spin indices are suppressed. 

The gauge transformation of the link variables is defined by 

Un,, - = glUn^n+fi- (2-9) 

where g n are elements of gauge group SU(iV) and t/ 9 denote the gauge transformation of U. 
By this gauge transformation, ip%(n, U) transform as 

V>f(n, 17) - ^(n, U°) = gl^(n, U), (2.10) 
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Then the vacuum overlap u(~\+)u is gauge invariant but the others are gauge dependent 
and are expressed as 



IJ9 



(±|±), = det ■ (2-li; 



The variation of the complex phase of the chiral determinant by gauge transformation 
gives the gauged Wess-Zumino effective action. 0) The corresponding quantity of the vacuum 
overlap formula is written by 

WHU,9) = f^^, (2.12) 

\ U g{±\±)l\ 

eM*r wz [U, 9 ]) = W - {UAYw+{u ^y (2-13) 

If all the link variables are set to be unity, it is expected that Eq. (2.13) gives the Wess- 
Zumino effective action (we call as Wess-Zumino term in the following). In this case the 
vacuum overlaps (|2.11|) are expressed by 

v(±\±)i = det |M(n,l)W,(n,l) , (2.14) 
and the Wess-Zumino term is defined by 

exp{% r wz \g\) = -— (2.15) 



For the continuum theory in two dimensions, it is well-known that the chiral determinant 
can be calculated exactly and the explicit expression can be obtained 0*. Next we verify 
that Eq. (2.15) has the correct continuum limit in two dimensions. 

2.2. Perturbative Calculation of Wess Zumino term 

First we examine ( |2.14| ) perturbatively. The group elements g n can be expressed as e l< ^ n 
by the variable (j) n which takes the value on the Lie algebra. We take this variables as the 
small quantum fluctuation from the unity (not some background) so that we can expand as 
follows: 

9n = e **» 1 + i<f> n + i(z n ) 2 + ■ ■ ■ . (2.16) 

Then, T wz \g\ can also be expanded in n and the continuum limit of the terms in the 
expansion series can be evaluated. 
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Let ^[g] denote the effective actions for Eq. ([2.14J) . Expanding by (f> n , we obtain the 
following result, 

'_pM-l 



^[9] = E I I tr(5g(k 1 )---6g(k M ))5o !kl+ ... +kM 

M= l 1V1 J ki,---,k M 

xC± M) {k ir --,k M ^), (2.17) 

Jp 

= jf tr (s ± (p)S ± (p -K x )--- S ± (p - K M -x)) , (2.18) 



where 



d 2 k 



(2.19) 



U (2tt) 2 ' 

Kj = E fei, ^0 = 0, (2.20) 

i=i 

(Wo" = (9nh - 5(*0 = E e~ ifc ' n ^„, (2.21) 

n 

0± = [1$ni>}], (2.22) 
^(Aj5 Q/3 + if ± (p) Q/3 ), (2.23) 



2 K 



H ± (p) is the momentum representation of Hamiltonians Eq. (2.3) in the free case and ijjp 
and Xp are the eigenfunctions and the eigenvalues of them. The explicit expressions are 
given in the appendix A. Then the Wess-Zumino term on the lattice is given by 



r wz [g] = lm{<P + [g]-<p-[g]) . (2.24) 



The oder contribution stems from M = 1 term, but Ct. and C ( ~ are equal to each other 
and cancel. The oder 2 contribution also vanishes since 2 term in M = 2 terms are real. 
Therefore the leading term is the order 3 and is given by the following, 



r wz [gf ] =i\l tr ($(kS(kS(k 3 )) So M 



O Jkl,k2,k3 

xlm (c+ (fci, fc 2 ) - C ( 7)(^i) • (2.25) 

where the explicit expressions of are also given in the Appendix A. In order to take the 
continuum limit, we introduce physical momenta Pi and express the dimensionless momenta 
ki as ki = Pid with the lattice spacing a. Then we expands the coefficient CtL in a and obtain 

r wz [gf ] = * 3 y pip2p3 7^7^ (^ tr (?(Pi)?(P»)fe)) 
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d 2 k 



8A 



±3 



(Cf(k)Cf{k)-Cf(k) 2 ) (B(k) T m ) 



Gvr 2 



+2 (C,{k)Ci v {k) - C u (k)C>»(k)) B{k)>* 



[2.26) 



[2.27) 



with />" = 9 M / . 

For the Wilson fermion, C and B are given by 



C M (fc) = sin/c^, = (1 - cosfc M ) . 



'2.28) 



Then the integrands of J 1 * 1 are complex expressions of sine and cosine and include two 
parameters, mo and r. It dose not seem to be able to accomplish the integrals analytically. 
By numerical calculation we obtained J + = it, J~ = for the region of m^jr , < m /r < 2. 
Therefore, when a goes to zero Eq. (2.26) becomes as follows: 



r wz [g\ 



(3) 



d 2 x e M „tr (<j)(x)d fl (j)(x)d u 4>(x)) . 



(2.29) 



a^O Y2-K 

This is identical with the leading term of the continuum Wess-Zumino term, 

r wz [g] c = d 3 y e ABC tr {g~ l d A g g- l d B g g- l d c g) (2.30) 

in the 0-expansion. 

2.3. Wess-Zumino- Witten model on the lattice 

Next we formulate a model with the Wess-Zumino term defined in the previous subsec- 
tions. The partition function of the model is defined by 



Z = I Vge- S[9] e inrwz[9] , 



(2.32) 



where S[g] is the action of SU(N) spin model: 
%] = ^£tr (^(n)^(n)t 



T2 J2 tr (p{n)g{n + fi) ] + g{n) ] g(n + //)) + const, (2.33) 



*) For the continuum vacuum overlap formula, 0) the same calculation can be performed. We obtain the 
similar expansion as Eq. (2.27) in the (£;/m)-expansion where m is Pauli-Villars mass. 



d 2 p 



(2.31) 



A P = [p 2 



It is easy to carry out this integral analytically and obtain the results J + = — J~ = ir/2. It also reproduce the 
correct Wess-Zumino term. The fact that J~ = ip lattice regularization is consistent with the calculation 
of the Chern-Simons current by Golterman et al.Q 
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and n is an arbitrary integer. This model is considered as the lattice counterpart of the 
two dimensional Wess-Zumino-Witten model.© In order to verify that the theory have the 
correct continuum limit, we perform the calculation of the Callan-Symanzik (3 function and 
show it is identical to that of the continuum model. 

To evaluate the (3 function, one may use the background field method, in which the field 
variables are separated into two parts such as 

5« = 9one Aw ", (2.34) 

with g 0n a smooth background and 7r n as small fluctuation around g 0n . Then one calculates 
the correction to the coefficient of the functional of go n , 

S c i[go] = -J2 tT (% Agon) 2 • (2.35) 

n,fi 

However, in our model, it is not easy to calculate the Wess-Zumino functional defined by 
Eq. (2.13) keeping the nontrivial background go n . So we first expand the total action in (f) n 
by which g n is expressed as e l ^ n . Then we shift <f> n around a background 0on n °t linearly 
but nonlinearly by applying the Hausdorff 's formula to Eq. (2.34) as follows: 

4>n = 4>0n + A7T n + -X[(f> 0n , ^n] 



rx (A[0on, [0o n, TTn] ] + A 2 [7r(n), [vr(n), <p 0n ] fj H . (2.36) 



12 

In this case, the lowest order of the functional (2.35) is 



Sd[4>o] = 51 tr (p^Ond^On) , (2.37) 



n,fi 



and it is sufficient for us to estimate the corrections to its coefficient. The one- loop order 
corresponds to 0(</>q, A 2 ). To obtain the propagator and vertices necessary to calculate the 
correction to this order, we first expand the actions S[g] and F wz \g\ in n as follow, 



A n,a LZ 

-T^(^n) 4 + ••■), (2-38) 
r wz [g] = I ai (^i)0 a2 (^)0 a3 (A;3)(5o, fcl+fe2+fc 3 

Jfci,fc 2 ,A;3 

xV® a ^ a3 (k 1 ,k 2 ) + ---, (2.39) 
V (3)a im{ku h) = _I/«««.] m ( C +(fc 1) k 2 ) - C-(k u k 2 )). (2.40) 
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Next we shift 4> n around 0o« as mentioned above and obtain the followings. 

S [(j>] = ^S cl [d>], (2-41) 
1 



where 



and 



SM = -J k <P a (-k)r(k)A(k)-\ (2.42) 
A{k)- l = £(l-cos^), (2.43) 



Si= I 0S 1 (^i)^ a2 (^)vf a3 (^)Vi+fe 2 +fe3 

xVl a ^' a3 (k 3 ), (2.44) 

5 2 = f C(^i)C(^)^ a2 (A;2)vf a4 (A;4)5o, fcl+ te + fe + fc4 

Jki,k2,k3,k4 

xV 2 aia2a - Aa4 {k u k 2 ), (2.45) 

5 3 = [ C(^i)C(^)^ 3 (A;3)7f a4 (A; 4 )5 0ifcl+fe2+fc 3 +fc4 

J k\ ,k2,ks,k4 

x(V 3 aia2asa4 (k 3 , k 4 ) + V 3 a3a4aia2 (k u k 2 ) + 4\/ 3 a3a2aia4 (A; 1 , fc 4 )), (2.46) 
54= / C(^i)C(^)^ a3 (A;3)^ 4 (A;4)5o, fcl+fe2+fc 3+fc4 

J fci ,k2,kz,ki 

xiAV^^ih, k 2 , fc 3 , fc 4 ) + 2y 4 a2aia3a4 (A; 1 , fc 2 , fc 3 , fc 4 )), (2.47) 

= 3A 2 / C(^i)^ a2 (^)^ a3 (A;3)5o, fcl+fe2+fc 3 

x^™ 3 ^,*;,), (2.48) 
r_, 2 = -3A 2 / C(^i)^ a2 (^)C(A;3)^ 3 (A;3)5o, fcl+fc2+fc 3 +fc4 

Jki,k2,kz 

xV£ )aiaae (k u k 2 )f a:i ' aie . (2.49) 

The vertex functions V's are listed in appendix B. Then the partition function is calculated 
as follows: 

^g-Sl- S2-S3-S4+in,r wz ,i+inr wZt 2 ^2 50) 

= Sofa] + ^((5i) 2 )o - (S 2 ) - (S 3 )o - (S 4 ) 

{X2,,) \(r^)\-inX^ni A )o, (2.51) 



2! 



where 



(O) = | P7rO(7r)e- 5 2°o 2) W. ( 2 .52) 
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The contribution of the original spin model is calculated as© 



-Z\(O)S d [0o]^ o - lna 2 m 2 ) S d [<f> ], (2.53) 



4 

where Ca is the seconde Casimir coefficient in the adjoint representation and m is the infrared 
cut-off mass. On the other hand, the contribution from the Wess-Zumino term is calculated 
as follows. It is easy to see the second term vanishes. The first term is given by 

- ^^((^' 2) ) 2 )o = \J r (-k)r (k) (h(k) + i 2 (k)) , (2.54) 



h(k) = -(^j C A J p A(p)A(k-p)C(k,-p) 2 , (2.55) 

h{k)=(^j C A J p A(p)A(k-p)C(k,-p)C(-k,k-p), (2.56) 

C(k,p) = lm (C+(k,p)-C+(k,p)) . (2.57) 

Toward the continuum limit we also scale the dimensionless external momenta k by a 
and consider the momenta ka near zero. Following Karsten and SmfiB, we can estimate 
the loop integral. The divergent contribution stems from the momenta near zero, so that 
we can use the continuum limit of the vertex C(k,p) obtained in the previous section. It is 
calculated as 



h (k) ~ I 2 {k) ~ JL ^-JcA^lna^. (2.58) 



Therefore the contribution from Wess-Zumino term is 



-^<(«„ = , (£) % (-L W ) S M + const. (2.59) 
Adding two parts, we obtain 

Z^exp^—SM + '-^J, (2.60) 



where A|, denotes the renormalized coupling constant given by 

1 1 / fn\ 2 \ 2 \ C A 1 , . 

rr '- d is" 1 (2 - 61) 



From this, (3 function is obtained as 
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For n > 0, the (3 function has an IR-fixed point at A 2 = — as expected (the factor 2 differ 
from Witten's results® because of the difference of the convention). For A 2 < — , this model 
is asymptotically free, so that it has a renormalization group invariant scale parameter, A 
given by 



A 




(2.63) 
(2.64) 



§3. Numerical Calculation of the Wess-Zumino term 

In this section, we show our results of the numerical calculation of the Wess-Zumino term. 
We consider the SU(2) chiral field and adopt the lattice action with the nearest neighbor 
coupling. We generate SU(2) chiral fields by the cluster algorithm.© Lattice size is set 
to 16 x 16. For most of the values of (3 = A examined, 500 configurations are generated. 
For (3 = 1.1, 1.5, we generate 4000 more configurations. With these configurations, the 
vacuum overlaps are calculated numerically through the LU decomposition. We measure the 
distribution of the Wess-Zumino term by dividing the range of the value [— it, ir] into bins 
with the width 0.1. We also calculate the standard deviation awz of the distribution of the 
Wess-Zumino term for each j3. 

First of all, we show the behavior of the Wess-Zumino term for randomly generated 
configurations of SU(2) chiral field in Figure 1. {(3 is set to 0.025.) The Wess-Zumino term 
distributes almost uniformly in the range [— tt, 7r]. In Figure 2, we show a typical configuration 
of the random SU(2) chiral field in momentum space. We can find that almost all modes of 
momenta contribute equally. The species doublers are then suspected to contribute to the 
complex phase. 

Next we show the distributions of the Wess-Zumino term for (3 = 0.5, 1.1 and 1.5 in 
Figure |[ We can clearly observe that for (3 = 1.1, the evaluated values of the Wess-Zumino 
term start to concentrate around zero and the distribution becomes like Gaussian. When j3 
increases, the distribution becomes sharp Gaussian. (3 dependence of the standard deviation 
&wz of the Wess-Zumino term distribution is plotted in Figure |j. The onset of the Gaussian 
distribution is around (3 = 0.9. Crossover is sharp from the strong coupling region where the 
Wess-Zumino term fluctuates hard. 

In Figure 5, we show a typical configuration of the SU(2) chiral field at /3 — 1.5 in 
momentum space. We can find that the modes with small momenta dominates. In this 
region of f3, the species doublers are almost suppressed. 
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Fig. 1. The distribution of the two-dimensional Wess-Zumino term for randomly generated SU(2) 
chiral fields. The range of the value [— 7r, tt] is divided into bins by the width 0.1. The normalized 
number of configurations for each bin is shown. L=16. 




Fig. 2. A typical configuration of random SU(2) chiral field in momentum space: the absolute 
value of the trace part is shown. L=16. 
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Fig. 3. The distribution of the two-dimensional Wess-Zumino term for j3 = 0.5, 1.1 and 1.5. 
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Fig. 4. The standard deviation <jwz of the distribution of the Wess-Zumino term for various j3. 




Fig. 5. A typical configuration of SU(2) chiral field at (3 = 1.5 in momentum space: the absolute 
value of the trace part is shown. L=16. 

In Figure || we show the susceptibility \ °f the SU(2) chiral model for L = 16 and 
L = 64. The scaling of the susceptibility starts around (3 = 0.9. This coincides with the 
onset of the Gaussian distribution of the Wess-Zumino term. 

In summary, the two-dimensional Wess-Zumino term defined by the vacuum overlap 
formula shows a sharp Gaussian distribution for the configurations of SU(2) chiral field in 
the scaling region. In this region, the configurations are smooth and the contribution of 
the species doublers is suppressed. In the strong coupling region, the Wess-Zumino term 
fluctuates hard and the species doublers' contribution is suspected to affect it. However the 
crossover is rather sharp from the strong coupling region to the weak coupling region. 

A few comments are in order for our numerical calculation. In Figure || we find that the 
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finite size effect on x is substantial for L = 16. We also suspect the finite size effect on the 
Wess-Zumino term. The onset of the Gaussian distribution of the Wess-Zumino term in the 
scaling region is also expected for larger lattice. However there can be much finite size effect 
in the values of <Jwz- We need further investigation on large lattice. 
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Fig. 6. Scaling of susceptibility of SU(2) chiral model. L=16 and L=64. 



§4. Numerical Evaluation of Observable in Wess-Zumino- Witten model 

In this section, we discuss the possibility of the numerical estimate of observables in the 
Wess-Zumino- Witten model by incorporating the Wess-Zumino action. The lattice Wess- 
Zumino- Witten model we are considering is defined with the complex action, 

Swzw[g] = S[g]- inr wz [g], (4.1) 

where S[g] and r wz [g] are given by Eq. (2.33) and by Eq. (2.15), respectively, n is the 
integer coupling constant of the Wess-Zumino term. In order to evaluate the observables in 
this model, we need to perform the following integral, 

(O) = J Vge- s[9]+inFwz[9] 0[g\ / Z, (4.2) 

Z = j Vge- s[9]+inFwz[9] . (4.3) 

In general, without the positivity of the measure, the Monte Carlo method cannot be 
applied. One possible way out is to incorporate the imaginary Wess-Zumino term into the 
observable and to perform the Monte Carlo integration only with the real part of the action. 
That is, we consider 

(0) = ((e mr ^O[g])) MC /z, (4.4) 
Z=((e mr ^)) MC , (4.5) 
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where (( ))mc denotes the Monte Carlo integration only with the real part of the action, 

{(X)) MC = fvge- s ^X[g). (4.6) 



Another possible method is to introduce the spectral functions© defined as follows: 

p{9) = Jvge~ s ^5(e-r wz [g}), (4.7) 

0(6) = J Vge- S ^5(6 - r wz [g])0(g), (4.8) 

(4.9) 

where the delta function is the periodic one. With these spectral functions, we can write the 
integrals Eq. (4.2) and Eq. (4.3) as 

(O) = r d6e ine O{6) lz, (4.10) 
Z = T d6e ine p(6). (4.11) 



The integration of 9 is then performed as the Riemann integral. 

Actually, the distribution of the Wess-Zumino term obtained in the previous section can 
be regarded as the spectral function, 

P (e) = ((5(6-r wz ))) MC , (4.i2) 

in the discrete approximation with the normalization, 

y £Ae P (e i ) = i. (4.i3) 

i 

We found that in the scaling region, it is given by a sharp Gaussian distribution with the 
width ay/z depending on /3, 

1 ~ i — (6»i+27rm) 2 

plfii) - -7= E e ^ • (4-14) 

VZrra W z 



m=—oo 



For such a Gaussian distribution of the imaginary part of the action, it is known that 
the first method mentioned above may work only for small n. © To see this, we perform the 
following test. We assume the Gaussian distribution for the Wess-Zumino term, Eq. Q4.14p . 



We take the square of the Wess-Zumino term as an observable. Since the spectral function 

*) It may be better to use the set methodil' for the precise evaluation of the spectral function. In this 
case, however, since the value of the Wess-Zumino term is in the compact region [— tt, tt], we can do without 
it if we have enough statistics. 
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is obtained by the Monte Carlo method, we simulate the situation by performing the Monte 
Carlo integration on 6. Thus we consider the following integrals for various n at fixed <Jwz, 



(4.15) 
(4.16) 



((r wz ) 2 ) = ((e m9 9 2 )) e /z, 
Z = ((e me ))e, 

where {{))$ denotes the Monte Carlo integration of 6 with the weight p(9), 

((X))o = d6p(6)X[d}. (4.17) 

J—n 

In Figure 7, we show the result. We found that for a W z = 0.5, up to n = 3 the numerical 
result reproduces the exact one within the acceptable error of three significant digits. For 
a wz = 0.2, n = 12. 

0.4 p , , , , , , , 

sigma_WZ = 0.2 » 

0.2 + % sigma_WZ = 0.5 + 




2 4 6 8 10 12 14 
Coupling Constant of Wess-Zumino term, n 

Fig. 7. Monte Carlo integration of 2 exp (in6) with the Gaussian weight. The width of the Gaus- 
sian weight awz is chosen as 0.2 and 0.5. The approximated analytical result for these awz, 
a wz(^ ~ a wz n2 )-> are a ^ so shown. 

On the other hand, it is also possible to perform the Riemann integral in Eq. (4.10) 
numerically for the same observable. We also found the upper bounds for n. In order that 
the numerical result reproduces the exact one within the acceptable error, we should have 
n < 12 for awz = 0.5 and n < 30 for a wz = 0.2. We have tried several schemes of higher 
order integration and obtained the similar result. 

Thus we obtained the upper bound for the accessible value of the coupling constant n of 
the Wess-Zumino term for given a wz . In other word, for a given n, there exists the upper 
bound of awz, with which the measurement of the observable can be performed by the direct 
Monte Carlo integration or by the Riemann integration of the spectral functions. In Figure 
9, we show the upper bounds of awz versus the coupling constant n of the Wess-Zumino 
term which is translated to the IR fixed point (5 C — From this figure, it seems possible 
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-12 



sigmaJ/VZ = 0.2 
sigma_WZ = 0.5 



5 10 15 20 25 30 

Coupling Constant of Wess-Zumino term, n 



Fig. 8. Riemann integration of 9 2 exp (inQ) with the Gaussian weight. The width of the Gaussian 
weight awz is chosen as 0.2 and 0.5. The approximated analytical result for these crwz, 
a wz(^ ~ a wz n2 )i are a ^ so shown. 



to simulate the model across the IR fixed point for f3 c > 1.3, (that is n > 17), if the spectral 
density method is used. On the contrary, by the direct Monte Carlo calculation, only the 
weak coupling region can be simulated. 



1.2 



0.8 



0.6 



0.4 



0.2 



Monte Carlo 
Riemann 
sigma_WZ 



0.6 



0.8 



1.2 1.4 
BETA 



1.6 



1.8 



Fig. 9. Upper bound on awz for given coupling constant n of the Wess-Zumino term, n is trans- 
lated to the value of the IR fixed point by the relation (3 C = The values of awz obtained 
in the previous section are also plotted. 

Our discussion given here is based on the numerical result shown in the previous section. 
Since the systematic and statistical errors are suspected in the result, the numbers of the 
upper bounds given in this discussion should not be taken literary. We believe, however, 
that this discussion gives a correct qualitative estimate of the possibility. 



5. Discussion and Conclusion 
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The vacuum overlap formula for the two-dimensional Wess-Zumino term was examined 
both perturbatively and nonperturbatively. We showed by the lattice perturbation theory 
that the formula correctly reproduces the Wess-Zumino term in the continuum limit and 
yields the IR fixed point in the beta function of the SU(2) chiral model. We calculated the 
vacuum overlaps numerically and found that the complex phase shows a sharp Gaussian 
distribution for the SU(2) chiral fields in the scaling region. In this region, the smooth 
configurations were obtained even in the symmetric phase due to asymptotic freedom. We 
also found that the crossover is rather sharp from the strong coupling region where the Wess- 
Zumino term fluctuates hard and the species doublers' contribution is suspected to affect 
it. 

It was also argued that, if we use the spectral density method, it seems possible to 
examine the region of the IR fixed point of the Wess-Zumino- Witten model numerically. 

Also from this study, we found that the asymptotically free coupling for gauge degree of 
freedom is able to reduce the gauge fluctuation. In two dimensions, the nearest neighbor 
coupling of the SU(2) chiral field plays such a role. In four dimensions, it is known in the 
continuum theory that the four-derivative coupling which is induced from the gauge fixing 
term is asymptotically free. It may be interesting to examine the gauge degree of freedom 
of the four- dimensional nonabelian gauge field with the gauge fixing term from this point of 
view. @) 

Acknowledgements 

The authors would like to thank H. Hata for enlightening discussions about the pertur- 
bative treatment of the non-linear sigma model. They are also grateful to W. Bock for his 
helpful advice about the Monte Carlo simulation with the Cluster algorithm. 

This work is supported by National Laboratory for High Energy Physics, as KEK Su- 
percomputer Project (Project No. 10). 

Appendix A 

Eigenvalues and Eigenfunctions of H ± (p) 

In this appendix we give the detail expression of the functions used in subsection 2.1. 
The momentum representation of Hamiltonians and the Wave functions are 

H ± (q)^(q)=X^(q), (A-l) 
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H%{P) 



B{p)Tm C(p) 
C(pY -B(p)±m 



Via 



1 / A± + B(p)Tm 



p 



— 



C*(p) 

C^p) 2 + (B(p) =f m ) 2 
2A± (Aj + J B(p) T m 



For Wilson fermion 



Ca»(p) = sin PM> 5 b) = r ^ (1 - cos k^) 



(A-2) 
(A-3) 
(A-4) 

(A-5) 
(A-6) 

(A-7) 



But one may take any regularization function as C(p) and -B(p). Imaginary part of C^, are 



IniCj^i,^) = / 
Jp 



B^(p; k u k 2 ) 

A\±\± \± 

A p-k 1 A p-k 1 -k 2 



B {3) (p; hi, k 2 ) = -£^u[Cn(p)C u (p - ki) (B(p -k x - k 2 ) =F m ) 

+ C„(p - fci)C v (p - fci - (B(p) T mo) 
+ C^p-k 1 -k 2 )C v (p)(B{p-k 1 )Tm G )}. 

Expanding on k±,k 2 , one obtain the expression of the integral J ± . 



(A-8) 



(A-9) 



Appendix B 

Vertices for One-Loop Calculation 



Here we give the vertex functions used in subsectection 2.3. 



^ 2 aia2a3a4 (fci, k 2 ) = f a ^ e f 



X 



1 -A(k 1 + k 2 )- 1 -^(A(k 1 )- 1 + A(k 2 )- 1 )y 



ya 1 a,2a 3 a 4 ^ \_ jaxa z e ja 2 a A e ^^(e*^ — X)(e ik ' M — 1) 



V, 



«1<12<13<14 



(hi, k 2 , k 3 , k 4 ) = tr (r^T a2 T as T ai ) 



x 



EC 



12 

l)(e 



i Am, 



l)(e 



l)(e 



1), 



VWwfa, k 2 ) = -l/™ Im (c+(h, k 2 ) - C+fa, k 2 )) . 



(B-l) 

(B-2) 
(B-3) 

(B-4) 
(B-5) 
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Appendix C 

Convention for SU(N) Lie Algebra 



ti(r a T b ) = l -5 ab , 



(C-l) 




(C-2) 
(C-3) 
(C-4) 



tr (T% d T b Ad ) = r d t cd = C A 5, 
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